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Abstract. We give Chebyshev-type quadrature formulas for certain new weight classes. 
These formulas are of highest possible degree when the number of nodes is a power of 2. 
We also describe the nodes in a constructive way, which is important for applications. One 
of our motivations to consider these type of problems is the Faraday cage phenomenon for 
discrete charges as discussed by J. Korevaar and his colleagues. 
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1. Introduction 

Let us denote by Pol{m) the family of polynomials of degree < m. 
Definition 1.1. Let p{x) be a nonnegative function on [-1, 1]. Denote by: 

the biggest natural number for which there exists a choice of distinct points {^/c})^^^^ c [-1, 1] 
so that the equality: 

(1.2) r P{x)p{x)dx = -fp{xu) 

is valid for an arbitrary polynomial P G Po/(M„(p)). 

Remark 1.3. A quadrature formula of the form (1.2) is called Chebyshev-type quadrature 
formula, and the points {xic}^^-^ c [-1, 1] are called the nodes of the quadrature formula. 
One can consider quadrature formulas of a more general form: 

P{x)p{x)dx = -yA,P{x,) 

where Ak are certain nonnegative numbers. These are known as Gauss type quadrature 
formulas. 

It is well known (see [10], p. 97) that a Gauss type quadrature formula with n nodes 
cannot be valid for an arbitrary P G Pol{2n). For completeness, we provide the proof of 
that fact here: 
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Remark 1.5. Let p{x) > be a nontrivial function on -1 < x < 1 . Then there are no points 
X]( G [-1, 1], k = 1, . . . ,n and numbers Aj^, k = 1,2, . . . ,n such that the formula (1.4) is valid 
for an arbitrary P G Pol{2n). 

Proof. Assume the converse, for some points xjt £ ^ = 1,2, ...,n and numbers 

A]^, k = 1,2, ... ,n the formula (1.4) is true. Then for the polynomial: 



k=l 



of degree 2n we must have: 



Po{x)p{x)dx = - Y AkPiXk) = 0. 
1 ^k~t 



□ 

Remark 1.6. Remark (1.5) implies thatM„(p) < 2n. 

The first Chebyshev-type formula was first discovered by Mehler in 1864 (see [10], p. 
185): 



Theorem 1.7 (Mehler). The formula 

is valid for an arbitrary polynomial P G Pol{2n - 1). 

On the other hand in 1875 K. Posse proved that (see [11] and [10], p.l86): 
Theorem 1.8. (Posse) If for a weight function 

p{x), X G [-1,1], 
we have that M„{p) = 2n - Ifor alln eN then: 

(1.9) p{x) ^ 



n 



The first weight function different from (1.9) for which M„{p) > n was given by UUman 
(see [5]): 

/ ^ 2 1 + bx ... 1 

where 

M„{[j) = n. 
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Later new examples of weights for whom M„(p) = n were found (see the survey article 
[6] and the references given there). These examples of weight functions are mainly of the 
following form: 

w{x) 
p{x) = , 

where w{x) is a certain positive analytic function on [-1, 1]. 

Examples for whom w has a singularity at have been found, too (see [4]). 

The main result of this article - theorem (4.2) provides Chebyshev-type quadrature 
formulas for certain new weight classes W„. 

All these classes include not only the weight: 

(1.10) p{x) ^ 



71 VT 



v2 



X^ 

(see remark (3.18)), but also many others (see remark (3.13)). 
Our result compares to the above-mentioned results by (3.19). 

The Chebyshev-type quadrature formulas for W„ are of highest possible degree when 
the number of nodes is a power of 2 (see remark (4.4)). Thus, even though (1.10) is the 
only weight for which the Chabyshev type quadrature formula is of the highest possible 
degree M„(p) = 2n - 1 (see theorem (1.8)), yet we prove that for every n which is a power 
of 2 there are many more weights for whom there exists a Chebyshev-type quadrature 
formula of the highest possible degree. 

Also note that the nodes in theorem (4.2) are described in a constructive way, which is 
important for applications. 



2. The Case p = 1/2 

Let us discuss the case p = 1/2 separately. S. Bernstein (see [10], page 197) proved that 
if n = 8 orn > 10 then one cannot choose points Xk,n, k = 1, 2, . . . , n in the segment [-1, 1] 
so that the following formula: 



If' iv- 

- P{x)dx = -)_^P{Xk,n) 

k=l 



is valid for all polynomials of degree n. 

S. Bernstein (see [9]) proved the following inequality M„(l/2): 

M„(l/2) < 7z V2n. 

A. Kuijlaars [3] (using the methods of S. Bernstein [8] ) proved that: 

Theorem 2.1. There exists an absolute constant c > such that: 

M„(l/2) >c^/n. 
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Today the interest into Chebyshev-type quadrature formulas is explained not only by 
numerous nontrivial conjectures relating to the issue (see [4],[1] for a list of open problems) 
but also a connection between the Faraday cage phenomenon for discrete charges and 
Chebyshev-type quadrature formulas, as explained by J. Korevaar and his colleagues (see 
[2]). For a more detailed explanation we refer the reader to the article [2]. 

This problem is also related to the discrepancy theory (see [7]). 



3. Auxiliary Constructions 

Let us denote by Xq = ~1 arid 



This sequence is increasing and 

lim x,j = 1. 

Definition 3.1. For an arbitrary natural n let us denote by Pni^) the polynomial of degree 
n, defined in the following way: 

1. we put: 

Po{x) = 1, P2{x) = 2x^ - 1, 

2. for an odd number q let: 

3. for a natural number n = 2^ let: 

4. For an arbitrary natural number n = 2^q, where q = l{mod2) and 3 < q let: 

Pnix) = {P2p{x)f . 

For n = 2^q, where q = \{mod2) denote: 

\Pn\ = V- 

Example 3.2. We have: 

Pq{x) = 1, Pi(x) = X, P2{x) = 2x^ - 1, Pjix) = x\ 

P^{x) = 2 (2x^ - if - 1, Psix) = x\ Peix) = (2x^ - if , 

P7ix)=x\ P8(x) = 2^2 (2x^-1)^ -ij -1 

Theorem 3.3. For an arbitrary natural n we have: 

max \Pk\ = [log, n]. 

0<k<n ^ 
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Proof. We present each number 1 < k < n n the form k = T'^qk, where = l{mod2). Then: 



max \Pk\ = max \pk\ = [log, n]. 

0<k<n 0<k<n ^ 



Example 3.4. We have: 

\P2^\=k, \P^,_,\=k-l. 

Theorem 3.5. For an arbitrary natural 1 <p we have: 

P2p{Xp) = -h P2p{Xp+l) = 0, P2P(1) = 1. 

On the interval [Xp, 1] the polynomial P2pix) is an increasing function. 
Definition 3.6. Let 

So : [xi, 1] — > [xo,Xi] , 

and for each x G [xi, 1], 

Sq(^x) = —X. 

For a natural n let us denote by: 

the one-to-one mapping, such that: 

Pi" (x) = -Pi" (S„(x)), x„+i < X < 1. 

We have: 

Sn{Xn+l) — Xji+l, Sh(1) = Xji, n = 1,2, . . . 



Example 3.7. We have: 



dSiiy) 

Si(y) = -y, = -1, X2<y< i. 



□ 



^ «y -1/2 

y + yi - y' ^53(y) y- yi-y' ^ 

Ssly) = ;= , — ; = , X4 < 1/ < 1- 

^ V2 rfy V2(i^ 

Definition 3.8. For an arbitrary integral number < n and a function /(x) defined on the 
interval x„ < x < 1, let us denote by Rnif) the function defined on the interval x„+i < y <1 
by: 

f (Sniy)) + fiy) 
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-1 -0,8 -0.6 -0.4 -0.2 0,2 0.4 0,6 0.8 



Figure 1. The mapping S2 

Remark 3.9. If the function / is constant on the interval [-1, 1] then for an arbitrary < n 
we have: 

Rn{f)iy) = fiyl on < y < 1. 

Remark 3.10. For any of the polynomials P)t(^)/ ^ = 0, 1,2, . . . which we constructed, we 
have: 

i^o(Pfc)(x) = P/c(x), k = Qimodl), 

and 

i^o(P^)(x) = 0, k = l{mod2). 

Remark 3.11. For any of the polynomials P2}'k{x), k,p = 1, 2, ... , which we constructed, we 
have: 

Rp{P2Pk){x) = PiPkix), k = Q{mod2), 

and 

Rp{Pvk){x) = 0, k = l{mod2). 

Definition 3.12. For natural n we'll denote by W,i the family of nonnegative functions 
p(x), -1 < X < 1 satisfying the following conditions: 

-1 

p{x)dx = 1 
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and 

dSkiy) 

piy) = ~P^^^^y'^^~j^' ^k<y<^, 

forfc = 0,l,2,...,n-l. 

Remark 3.13. Let us note that for an arbitrary weight function p{x) G W„ we have: 
(3.14) r p{y)dy = 2-\ 

Jx„ 

Moreover, an arbitrary nonnegative function defined on [x„, 1] satisfying the condition 
(3.14), coincides with a weight function from W„ on 1]. 

Example 3.15. Note, that: 

So(i/) = -y, xi<y< 1, 
therefore for the weight function p G Wi we get the condition: 

piy) = pi-yl xo<y<i. 

Example 3.16. We have: 

Siiy) = ^1 - y^> X2 < y < 1, 

therefore for the weight function p G W2 we get the conditions: 

piy) = pi-y)' ^0 = -1 < y < 1, 

and 

piy) = p[^'^-y^) ^jY—yi ' ^i<i/^i' 

In particular, the function: 

piy) = ^ 



satisfies these conditions. 
Example 3.17. We have: 

Siiy) = ^ ^ ^ , X3 < y < 1, 
V2 

therefore for the weight function p G W3 we get the conditions: 

piy) = pi-y)r xo < y < 1, 

and 

y 

, xi < y < 1, 



p(y) = p(V^-y')-^ 



y' 
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V2 



p(y) = p 

In particular, the following function: 

piy) 
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(y+ Vl^U- 



V2(i - y') 
1 



■, X2 < y <1. 



71 

satisfies all the three conditions. 

Remark 3.18. By a reasoning similar to (3.16) and (3.17) one can prove that: 

1 



piy) = 



e W„ for n = 1, 2, , 



71 

Remark 3.19. The class W„ coincides with the family of non-negative functions p{x) whose 
integral is equal to 1 and who permit the following representation: 

wix) 

where w{x) satisfies the conditions: 

w{x) = w(S)t(x)), for all = 0, 1, . . . , n - 1, and Xk+i < x < 1. 

4. Main Results 

Theorem 4.1. Lei n>lbea natural number. Let p{x) e W„. Then for an arbitrary function f{x) 
we have: 

r f{x)p{x)dx = r r {Rn-^o...oR,){f){x)-p{x)dx. 

J-1 Jx„ 

Proof If p{x) G Wi we have p{y) = p{-y) so, 

f{x)p{x)dx = I f{x)p{x)dx + I f{x)p{x)dx = 
1 J-i Jo 

= ^ /(^)p(^)rfx ~ X ^ " 

= r (/(x) + /(-x))p(x)^ix = 2 r j^o(/)(y)p(y)^^y. 

Jo Jxi 

If p{x) G W2 then we have p{y) = p{-y) and 

dSi{y) 

piy) = -p (Si(y)) ^2 < y < 1- 
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Thus, for an arbitrary function f{x) we have: 

f{x)p{x)dx = 2 j Ro{f){x)p{x)dx = 2 I Ro{f){x)p{x)dx + 2 I Ro{f){x)p{x)dx = 

1 Jxi Jxi Jx2 

= 2 J Ro{f){x)p{x)dx - 2 £ Roif) (Si(y)) p (S,{y)) ^^dy = 
= 2 f (Roif) (Si(y)) + Ro(/)(y)) p{y)dy = 2' f (Ri o Ro){f){y)p{y)dy. 

Jx2 J X2 

We prove the theorem for the other values of n in an analogous way. □ 
Theorem 4.2. Let p{x) G W„. Let X = {tji^l-^ be a certain set of points, for whom: 

-1 < t^ < t2 < ■ ■ ■ < tm < 1. 

Let X n [x„, 1] ^ and 

S)c (X n [X)c+i, 1]) = X n [X)t/^)c+i] / for all k = 0,1,2, ... ,n - 1. 
Then for an arbitrary polynomial P{x) e Pol{2" - 1) we have: 



P{x)pix)dx = - y P{tk). 
1 m ^ 

1 k=l 



Proof. Remark (3.11) implies that: 

(R„_io...oJ^o)P^ = 0, for fc = l,...,2"-l. 

So, for each polynomial P{x) e Pol (2" - 1) we have (Rn-i o ■■■ o Rq) P = const. 
Thus, by theorem (4.1) we have: 

r P{x)pix)dx = (Rn-i o • • • o Ro){P){t) ■ 2" f p{x)dx, 

J-l Jx„ 

for any x„ < t < 1. Further, by (3.13) we get: 

P{x)p{x)dx = (R,_i o . . . o Ro){Pm, 
for any x„ < i < 1. This is why, 

P{x)p{x)dx= iRn-lO---oRo){P){tj). 
^ I I n/ 1 I ;(.^^<t.<i 

At this moment let's consider the particular case n = 1. From (3.8) we have: 

1 „ . 1 



1 1 



|Xn[0,l]| ^ " 2|Xn[0,l]| L-i \ ^ >' ^ m 

' I- ' J I 0<fy<l ' I- ' J I o<ty<l k=l 
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For the other values of n the equality: 

^^^^^ 2:(R,.,o...„R„)(P)(,,)= 2^Pfe). 

' t «' J I xn<tj<\ k=l 

is proved in an analogous way. 

□ 

Remark 4.3. Theorem (4.2) implies that for any p G Wn we have: 

m = \X\> 2"-\ 

If we additionally impose that: 

then we get that \X\ = m = 2"~^ and the quadrature formula is valid for an arbitrary 
poynomial P{x) G Pol (2" - 1). Hence, 

M2,.-i(p) >2"-l 

or 

M„,(p)>2m-1, for m = 2""^ 
But taking into account the remark (1.6), we get that: 

M„,(p) = 2m - 1, for m = 2""\ 

Remark 4.4. Remark (4.3) shows that we found a Chebyshev - type quadrature formula for 
p G Wn of the highest possible degree 2" - 1, where the number of nodes is m = 2""^. 
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